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ABSTRACT 

Microlensing event OGLE-2002-BLG-055 is characterized by a smooth, sUghtly asymmetric single- 
lens curve with an isolated, secure data point that is ~ 0.6 magnitudes brighter than neighboring 
points separated by a few days. It was previously suggested that the single deviant data point and 
global asymmetry were best explained by a planetary companion to the primary lens with mass ratio 
q = 10"'^ — 10~^, and parallax effects induced by the motion of the Earth. We revisit the interpretation 
of OGLE-2002-BLG-055, and show that the data can be explained by wide variety of models. We 
find that the deviant data point can be fit by a large number of qualitatively different binary-lens 
models whose mass ratios range, at the ^ Scr level, from q ~ 10""' to ~ 10~^. This range is consistent 
with a planet, brown dwarf, or M-dwarf companion for reasonable primary masses of M > O.SMq. 
A subset of these binary-lens fits consist of a family of continuously degenerate models whose mass 
ratios differ by an order-of-magnitude, but whose light curves differ by < 2% for the majority of the 
perturbation. The deviant data point can also be explained by a binary companion to the source with 
secondary /primary flux ratio of ~ 1%. This model has the added appeal that the global asymmetry 
is naturally explained by the acceleration of the primary induced by the secondary. The binary- 
source model yields a measurement of the Einstein ring radius projected onto the source plane of 
fE = 1.87 ±0.40 AU. OGLE-2002-BLG-055 is an extreme example that illustrates the difficulties and 
degeneracies inherent in the interpretation of weakly perturbed and/or poorly sampled microlensing 
light curves. 

Subject headings: gravitational lensing ~ planetary systems 



1. INTRODUCTION 

Planetary companions to Galactic disk and bulge mi- 
crolens stars can be discovered via the short duration 
perturbation they create to the smooth light curve in- 
duced by the parent star l)Mao fc Paczvh ski 1991). The 
majority of these perturbations are expected to be rela- 
tively simple and grossly characterized by three observ- 
ables: the duration, peak time, and magnitude of the 
perturbation. In the ideal scenario, these three observ- 
ables are simply related to the three parameters describ- 
ing the planetary system: the planet/star mass ratio, 
the instantaneous projected separation in units of the 
angular Einstein ring radius, and the angle of the source 
traje ctory relative to the planet /star axis ( Gould & Locb 
Il992j) . Unfortunately, reality is a bit more complicated, 
and a number of degeneracies have been identified which 
can hamper the ability to in fer these parameters in prac- 
tice. iGaudi fc Gouldl ^^13) demonstrated that there ex- 
ists an ambiguity in the physical mechanism that sets the 
width of low-amplitude perturbations which can result in 
an order-of-magnitude uncertainty in the inferred mass 
ratio. Gaudi (1998) pointed out that a subset of binary 
sources with extreme flux ratios can reproduce the dura- 
tion and magnitude of a subse t of pla netary microlensing 
perturbations, although iHarJ l|2002f> demonstrated that 
this degeneracy could be resolved w ith astrometric ob- 
servations during the perturbation. iGriest fc Safizadehl 
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l)1998() discuss a two-fold discrete degeneracy in the pro- 
jected separation prevalent in high-magnification plane- 
tary events. Along with these anticipated degeneracies, a 
few have been uncovered in the process of det ailed mod- 
eling of observed events. iBennett et alJ l)1999fl invoked a 
planetary companion to explain a short-duration devia- 
tion seen on a clos e binary-lens l ight curve. However, it 
was later shown bv lAlbrow et al.l 1^200 0) that this pertur- 
bation could also be fit by one of the secondary caustics 
of the close binary lens, when rotation of the binary is 
considered. fCaudi et al. (2002) found a weakly asym- 
metric event that could be equally well-explained by a 
planetary companion, or parallax deviations arising from 
the motion of the Earth. Many of these degeneracies are 
'accidental,' in the sense that they arise from chance sim- 
ilarities between deviations caused by different physical 
situations, rather than by intrinsic degeneracies in the 
lens equation itself. They are therefore generally only 
approximate degeneracies, and can be resolved with ac- 
curate, well-sampled light curves. Given the short dura- 
tion and unpredictability of planetary deviations, dense, 
continuous and accurate light curve coverage is necessary 
to both detect and accurately characterize planetary mi- 
crolensing perturbations. 

A lensing star with a planetary companion is just an 
extreme limit of a binary-lens. As discussed by numerous 
authors, binary l enses are the mselve s subject to numer- 
ous degeneracies llAlbrow et al.»1999clJaroszvnski fc Mad 
120011: iMao fc Pi Stefanclll995H . some of which are roote d 
in symmetries in the lens equation itself fDominik"199 
and are therefore nearly perfect fAfonso et al. 200 
lAlb row et ani2002|) . Binary lenses in which the source 
does not cross any caustics can also be confused with bi- 
ktKEiH^'I sources. This is especially problematic when only 
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single-band photometry is available. This may partially 
account for the fact that, although they were predicted 
to be plentiful l|Griest fc Hulll99^ . only one candidate 
binary-source lensing event has been identifie d^, event 
OGLE-2003-BLG-095 ( Collinge"1994) . Indeed, 'CoUingd 
(fl994) found that the binary-source model for OGLE- 
2003-BLG-095 is only preferred over a binary-lens model 
at the ~ 3(7 level. 

Source and lens binarity are not the only regimes 
where degeneracies are plentiful; global deviations from 
the fid ucial point-source, point-lens, uniform motion (i.e. 
iPaczvh ski 198^ light curve have also been found to 
be subject to degeneracies. Such deviations come in 
a variety of forms. The motion of the Earth pro- 
duces departures from uniform relative motion which 
can induce observable deviations from the standard 
lightcurve form. These deviations can be quite dra- 
matic for even ts with timescales <e of order or larger 
than a year ijSmith et al.l I2002|) . However, in the 
more usual case where ^ yr, the effect of the 
motion of the Earth can be approximated by a con- 
stant acceleration, which results in deviations that can 
be symmetric or asymm etric w ith respect to the peak 
of the event iGould. M iralda- Escude. fc Bahcalil Il994t 
iSmith. Mao, fc Paczvhskii i2003jl . although asymmetric 
deviations are generally easi er to recognize. Unfortu- 
nately , as demonstrated by ISmith. Mao, fc Paczvhskii 
l)2003j) . any such weak parallax deviation can also be 
explained by acceleration of the source, thus mak- 
ing the parallax interpretation non-unique for such 
short timescale events. Smith. Mao, & Paczvhski (2003) 
furthermore demonstrated that constant-acceleration 
events are subject to a two-fold discrete degeneracy be- 
tween the magnitude an d direction of the acceleration 
and the event timescale. iGouldl l)2003j) showed that, in 
fact, light curve degeneracies extend to even higher or- 
der, and are present when one takes into account not only 
the acceleration, but also the jerk. This jerk-parallax 
degeneracy has b een seen in one event toward the bulge 
ijPark et al J 120041) . and has been invoked to resolve the 
discrepancy between the photometric and microlensing 
mass determinations of the microlens in the Large Mag- 
ellanic Cloud event MACHO-LMC-5 (|Goi]ld .2003V 

Galactic bulge microlensing event OGLE-2002-BLG- 
055 exhibits both a global asymmetry and single deviant 
data point, fjaroszvhski & Paczvhski (2002) first consid- 
ered the interpretation of this event; they fit binary- lens 
models and included non-uniform motion caused by ei- 
ther parallax or arbitrary uniform acceleration. They 
argued that, although both uniform acceleration and 
parallax explain the global deviations equally well, the 
event was most naturally explained by parallax-induced 
deviations, and that the single deviant point was best 
fit by a binary- lens with mass ratio 10^'^ — 10^^. This 
mass ratio implies a companion in the Jupiter mass range 
for likely primary masses, and they therefore concluded 
that OGLE-2002-BLG-055 was a possible planet candi- 
date. Here we revisit the interpretation of this event. 
We demonstrate that, in fact, there are many possible 
interpretations of OGLE-2002-BLG-055. We find that 
the short-timescale perturbation could arise from a stel- 

* See IDominiM IT993) and IHan fc Jeonel IT^M] for additional 
discussions of the apparent lack of binary-source events. 



15.5 



0) 



III 

i 


1 1 

15.4 




1 






1 1 








15.6 






- 


- A) 


15.8 






_ 


/ \ 
J 


16 


-V / 








\ 16.2 












" , 1 








J 2420 


2440 












1 1 1 1 


1 1 


1 


'i 





2400 2500 2600 

HJD-2450000. 

Fig. 1. — The points with errorbars show the Hght curve of 
OGLE-2002-BLG-055. The long-dashed magenta line shows the 
best single-lens, single-source, constant-velocity model fit to the 
data, with the single high point and two neighboring points re- 
moved. The short-dashed blue line shows the best single-lens, 
single-source, constant-acceleration model fit to the same dataset. 
The dotted green line shows the best single-lens, binary-source, 
constant-acceleration model to the entire data set. The solid red 
line shows the best single-source, binary-lens, constant-acceleration 
model to the entire dataset. The inset shows a close-up near the 
deviation. 



lar companion to the source, or a stellar, brown-dwarf, or 
planetary companion to the lens. The global asymmetry 
could arise from parallax deviations, or from accelera- 
tion of the source induced by its companion. All of these 
interpretations are indistinguishable at the ~ 3cr level. 
Therefore, the correct interpretation of this event is un- 
clear, and it serves as a particularly extreme reminder 
of the degeneracies involved with poorly sampled and 
weakly-perturbed microlensing light curves. 

2. DATA 

The Galactic bulge microlensing event OGLE-2002- 
BLG-055 was observed a s part of the third phase of 
the OGLE collaboration l|Udalsk i 20031^, and alerted 
during the 2002 bulge season with the Early Warn- 
ing System (EWS)^ The OGLE data consist of 106 I- 
band data points, with 9 points taken during the 2001 
season, 47 during the 2002 season, and 50 taken dur- 
ing the 2003 season. The data covering the primary 
event are shown in Figure ^ All errors are scaled by 
a factor of 1.31, as determined from the 2003 baseline 
data. This scaling may be somewhat underestimated, 
as the best models found here have unreasonably large 
values of x^, although they appear to faithfully repro- 
duce all the features of the event. The majority of 
the data for OGLE-2002-BLG-055 appear to follow the 
usual single-lens form, with the obvious exception of 
the data point at HJD' = HJD - 2450000. = 2424.9, 
which has / — 15.449 ± 0.009, highly discrepant from 

® http:/ /www. astrouw.edu.pl/~ftp/ogle/ogle3/ews/ews. html 
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the n eiRhborinK points at / ^ 16 separated by 4 — 8 
days. I.Taroszvnski fc Paczvn ski (2003 report that this 
data point is secure, as evidenced by direct inspection of 
the raw image. 

Unfortunately, there is only /-band data available for 
OGLE-2002-BLG-055, which hinders the interpretation 
of the model fits to the data for several reasons. First, 
without color information during the event, it is impos- 
sible to determine separately the color of the source and 
blend. Thus it is not possible to compare the derived 
source and blend colors and magnitudes to a local color- 
magnitude diagram, which provides an important test 
of the viability of a given model. In fact, as this field 
was not observed during the second phase of the OGLE 
collaboration, there is no TZ-band data available for the 
source field at all. Thus it is also not possible to deter- 
mine the local /-band e xtinction At via the usual method 
of 'cl ump-calibration' ijStaneU [19961: iWozniak fc Stanekl 
Il996l^ . There are several publicly-availab le ext i nction 
maps for the Galactic bulge. The maps of ' Sumil l|20fl4f) 
are based on OGLE 11 data, and thus do not cover the 
OGLE -2002-BLG-055 field. iPopowski. Cook, fc Beckal 
l)2003t) have derived extinction maps based on MACHO 
collaboration data. We estimate the extinction at the 
location of the source (R.A. = 17''59'"40f93, Dec. = 
-27°07'18."2(J2000); I = 3.14115,6 = -1.75827) as the 
mean of the t wo closest points (at distanc es of 0.12° and 
0.16°) in the iPopowski. Cook, fc Beckeil l)2003|) extinc- 
tion maps, weighted by the inverse squared distance. 
This yields Ay = 3.20, or A j = 1.57, where we have 
adopted Ai/Av = 0.49 ( Sui^ l200l FUdalski 2003a!). 

3. MODELING OGLE-2002-BLG-055: TECHNICAL DETAILS 

We fit the flux as a function of time to the standard 
form. 

Fit) - F,A{t) + Ffc, (1) 

where Fs is the flux of the source, and Ff, is the flux of any 
unrelated blended light. Here A{t) is the magnification 
as a function of time t. 

For a single lens (SL) and single source (SS), the mag- 
nification is, 



Ao[u{t)] = 



u\Jx? 



(2) 



where u(t) is the angular separation between the lens and 
source in units of the angular Einstein ring radius, 

= \J nM-K^ch 1^ = n .-TT - (3) 
c^AU 

Here tt^cI is the lens-source relative parallax, and M is 
the mass of the lens. 

The angular separation between the source and lens 
can be written as, 

u'{t)^x\t)+y^{t). (4) 

We consider two different models for the relative lens- 
source motion. For constant relative velocity (CV) be- 
tween the source, lens, and observer, 

x{t) = r = ^ yit) = uo (5) 

where = O^/fi^ci is the Einstein timescale, fi^ci is the 
relative lens-source proper motion, is the time of clos- 
est approach between lens and source, and uq is the im- 
pact parameter in units of ^e- For constant acceleration 
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Fig. 2. — In each panel, the points with errorbars show a close-up 
of the light curve of OGLE-2002-BLG-055 for ~ 15 days before and 
after the deviant data point. The solid curve in each panel shows a 
different single-source, binary-lens, constant-acceleration fit. The 
mass ratio and of each fit are also indicated. The dotted curve 
shows the best single-source, single-lens, constant-acceleration fit 
with three data points removed, and is the same as the dashed 
curve in Figure The bottom two panels show a continuum of 
fits, with the range of mass ratios and values indicated. 



(CA), the two components of the angu lar separation be- 
come l|Smith. Mao, fc Paczvhskill2003D . 



x{t) 



y{t) = Uo + 2 A'^^ 



(6) 



Here Ax and Ay are the two components of the relative 
source- lens angular acceleration in units of Note 
that and /iroi are now defined at time t^. 

For a SL and binary source (BS), the magnification can 
be written 

^BS [t) - ^0 [ui (t)] + e^o [U2 it)] , (7) 

where e is the flux ratio between the secondary and pri- 
mary source. Here ui and U2 are the angular separations 
between the primary and secondary source and the lens, 
respectively. 
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For a bound binary-source, we expand the vector an- 
gular separation of the primary and secondary around 
the time of maximum magnification of the primary to, 
and keep terms up to the acceleration. This yields ui{t) 
given by Eq. with to and uq the time of maximum 
magnification and impact parameter of the primary, and 

/ I \ 2 , ; / < \ ^ •^V 2 , ; 

X2(t) = Ti - - — Ti + bx, y2(t) = U0-- Ti + by, 

where ti — {t — tQ)/t^, and is the mass ratio of the 
binary-source system, and b = {bx,by) is the vector an- 
gular separation from the primary to the secondary in 
units of 0E, which has the same direction as the accel- 
eration of the primary. Unfortunately, this parameteri- 
zation is not ideal for the current event because bx and 
by are not directly constrained by the data, but rather 
the magnification of the secondary source at the time 
of the deviant data point, A{td) — e/u2{td)- Therefore, 
adopting Eq. |S1 would lead to correlations between the 
parameters Ax, Ay, to, uq, bx, and by, which would lead 
to complications with fitting and error determination. 
Instead, we adopt X2{t) = {t — to,2)/iE and 2/2(0 = mo,2- 
This assumes that the acceleration of the secondary is 
negligible, and that the direction of the acceleration of 
the primary is not necessarily the same as b. The ac- 
celeration of the secondary has two effects. First, the 
effective timescale at the time of the perturbation may 
be different from ^e- Second, the deviation from uni- 
form velocity induces a deviation in the light curve. The 
former effect is essentially unobservable, because the dif- 
ference in can be absorbed into e and uo,2 ijGaudil 
119981 IcTouldl IT996) . For typical parameters, the devia- 
tion in the light curve due to acceleration is < 1% while 
the magnification of the secondary is significant. This is 
of order or smaller than the errors, and thus is negligible. 
The alignment between the source separation vector and 
acceleration vector as determined from the fit provides 
an important test of the binary-source interpretation. 

For a binary lens (BL), the magnification is specified 
by three parameters: the mass ratio q of the lens, and 
instantaneous angular separation d between the two com- 
ponents in units of and the angle a between the vector 
pointing from primary to secondary and the direction of 
/Xrei. Here we adopt a coordinate system such that the 
origin of the binary- lens is at a distance of qd~^{l + q)~^ 
from the primary along the axis toward the secondary. 
For BL models with small mass ratios, it is also necessary 
to consider the effect the finite size of the source star 6** 
on the magnification, parameterized by = 9^/9e- 

We perform all fits in flux rather than magnitude space, 
thus allowing the analytic determination of the values of 
Fg and Fj, for a given model A{t) that minimize the usual 
goodness-of-fit statistic. All other parameters are 
minimized using the downhill-simplex rout i ne AM OEBA 
(jPress. Teukolskv. Vetterling. fc Flannervl Il992|) . using 
many different initial trial values of the parameters as 
seeds. This procedure works well for all models except 
the BL model, which generally produces poorly-behaved 
surfaces that are not well explored by downhill- 
simplex routines. 

In order to approximately minimize BL models, and 
to calculate errors on fit parameters for other models, we 



use a Markov Chain Monte Carlo (MCMC)*^. At each 
step in the chain, we vary a single fit parameter by 
drawing a random Gaussian deviate with zero mean, 
and dispersion set to appropriately sample the likeli- 
hood surface. We then calculate the relative likelihood 
Crei — exp(— between the new step, and the old 
position. If Crei > 1, the step is accepted. If £rei < 1, 
then we draw a random uniform deviate between and 
1. If this deviate is < Crei, then the new step is accepted, 
otherwise it is rejected. When generated in this way, the 
resulting distribution of parameters in the Markov Chain 
is proportional to the posterior probability distribution, 
provided that the chain has converged. Although we do 
not rigorously test for convergence, we confirm that the 
chains have been run for a sufficient number of steps, by 
visually inspecting different chains started at different 
points in parameter space, and ensuring that these have 
properly mixed. 

We define the la errors on a given parameter as the 
projection of the Ax^ = 1 contour onto that param- 
eter axis. We do not attempt to calculate errors for 
the BL models, primarily because we find that there are 
large, disconnected regions of parameter space that con- 
tain continua of degenerate models. Defining confidence 
regions in the presence of such degeneracies is difficult 
and not very meaningful. 

4. MODELING OGLE-2002-BLG-055: RESULTS 

Table 1 shows fit parameters and, where appropriate, 
Icr errors for the four different classes of models consid- 
ered here. The models are labeled by the lens multiplic- 
ity (SL=single lens, BL=binary lens), source multiplic- 
ity (SS=single source, BS=binary source), and source 
motion (CV=constant velocity, CA=constant accelera- 
tion). Also indicated is the dataset used; no number 
refers to the entire dataset, '1' indicates the dataset 
with the single deviant data point at td = 2424.9 re- 
moved, '3' refers the dataset with points at t — 2416.81, 
2424.9, and 2428.89 removed. For the constant accelera- 
tion fitSj_tlier£_a£e_genei^^ two degenerate mod- 
els (jSmith. Mao, fc Paczvhskil200,H) . these are labeled 'p' 
and 'n' for the fits with positive and negative Ay, respec- 
tively. The binary-lens fits are labeled as (a,b,c,d,e,f), 
which corre spond to the panel labels i n Figu re El 

Following iJaroszvhski fc Paczvnskil l)2002f) , we first fit 
OGLE-2002-BLG-055 to a standard single-lens, single- 
source, constant velocity model after removing the sin- 
gle deviant data point at td — 2424.9. This model is 
clearly a poor fit to the data, yielding — 180 for 100 
degrees of freedom (dof). Removing two additional data 
points on either side of the deviant data point improves 
the fit by Ax^ — 12 for two fewer dof. This implies 
that the two data points on either side of the point at td 
likely deviate significantly from the smooth underlying 
model, supporting the reality of the single data point, 
and pointing toward a perturbation that has a timescale 
that is of order the interval between the td and neigh- 
boring data points, ^ 4 days. All constant-velocity fits 
are quite poor, deviating significantly from the data dur- 
ing the rising part of the primary event (see Figure^, 
likely indicating a significant acceleration of the observer, 

^ See I Verde et alJI2003 for a particularly clear and concise ex- 
planation of using MCMC in a practical setting. 
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source, or lens. 

Including a constant acceleration improves the fits con- 
siderably. For the data set with one data point re- 
moved, we find — 149 for 98 dof. When three 
data points are removed, the fit improves by Ax^ = 
15.7 for 2 less degrees of freedom. As anticipated by 
ISmith. Mao, fc Paczvhskil l|2003t) , we find two degenerate 
fits for each constant acceleration model, with different 
values of ^x, ^md Ay, but identical values of the other 
parameters and x^- Including the constant acceleration 
terms increases the inferred timescale considerably. The 
value of for the best-fit constant acceleration model 
with three data points removed is unacceptably high, 
despite the fact that the model appears to reproduce 
the primary features of dataset quite well (see Figure 
^ . Inspection of the light curve during the microlensing 
event, while the source is significantly magnified, reveals 
short timescale (< 1 day) scatter, with amplitude that 
is smaller than the typical size of the photometric errors 
at baseline. This scatter may be intrinsic to the source. 
Regardless of the cause, it appears that the error scaling 
derived from the baseline points is probably underesti- 
mated. 

A binary-source model can produce a large, short- 
timescale deviation like that seen in the dataset of 
OGLE-2002-BLG-055, provided that e < 1 and uo,2 < 
1. Figure El shows the best-fit binary-source models, 
which have — 142 for 96 degrees-of-freedom. This 
is Ax^ — 11 larger than the single-source constant- 
acceleration fit with three data points removed. This 
additional arises from the inability of the binary- 
source model to simultaneously fit the two data points 
immediately following td- As with the single-source 
models, there are two binary-source models with es- 
sentially equal x^ due to the constant-acceleration de- 
generacy. As anticipated, the fits yield small values 
for the flux ratio, e ~ 1%. In this regime, the rela- 
tions between the binary-source parameters e, uq, and 
^0 and the salien t features of the perturbation are ex- 
ceedingly simple ( IGaudilll998(l . Note that the param- 
eters of the primary (to , i e , , Ax , Ay , Fg, Ff,), which 
are constrained by the data away from the perturba- 
tion, are essentially identical to the parameters found for 
the constant-acceleration fits with three data points re- 
moved. The duration of the perturbation is constrained 
by the points neighboring td- Combined with the param- 
eters of the primary, this duration yields e. The excess 
flux at time td is AF{td) — Fse/u2{td), and thus yields 
iJ'2itd), which is much better constrained than to, 2 and 
uq 2 separately. We find U2{td) = 0.0025 ± 0.0008 and 
U2{td) = 0.0018 ± 0.0007 for models SL.BS.CA.n and 
SL.BS.CA.p, respectively. Note that, since the relations 
between the perturbation observables and the binary- 
source parameters e, to,2, wo,2 depend on ^e, which differs 
between the two constant-acceleration fits, the binary- 
source parameters also differ between the two fits. 

Binary-lens models with extreme mass ratios q can also 
yield large, short-timescale deviations. In the limit of 
q 1, and under some simplifying assumptions, the 
relation between the gross features of the perturbation 
and the binary pa rameters q, d, a also t ake on a rela- 
tively simple form ijGaudi fc Gouldlll997t IGould fc Loebl 
Il992j) . The angle of the trajectory is roughly a ~ 



tan^"'^[iio/'''(^d)]- The binary separation is d ^ 0.5[ud + 
+ 4] , where Ud ~ V^i^d)^ + Wq- Finally, the 
mass ratio is q ^ ^jAt/t^. For fits SL.SS.CA.3n and 
SL.SS.CA.3p, these expressions yield a ~ 0.6 — 0.8, 
d ^ 1.1, and q ^ 10^"^. Of course, these formulas are 
only very rough approximations, but the implied param- 
eters are good starting guesses. We therefore start with 
these parameter ranges, and first alter the parameters by 
hand until we find a number of promising approximate 
fits. These fits are then used as seeds for an automated 
point-source binary-lens minimization routine. All viable 
point-source fits are then collected, and the parameters 
are used as starting guesses for a second minimization, 
using a finite source of angular size 

=6/.asl0-°-2(^«-i4-32). (9) 

Here /q is the dereddened magnitude of the source 
as derived from the point-source fits. Equation lO 
is derived from the color-surface brightness relation of 
Ivan Belli (|1999() . assuming that the {V — I)o color of 
the source is the same as the red clump. Standard 
models of the Galaxy ijHan fc Gouldl IToM 120031 pre- 
dict (e'E) = 240ti!^|^/ias (bulge), and (6'e) = 330tJ^|^Aias 
(disk). These values of 9b yield source sizes — 6'*/6'e 
that are too large to reproduce the data for OGLE-2002- 
BLG-055 with small mass ratios q. Although this implies 
that such fits are somewhat disfavored, it is not actually 
possible to rule out these models with such an argument. 
We therefore adopt a value of 9e = 550/ias. Approx- 
imately 6% of bulge lenses and 25% of disk lenses are 
expected to have values of 9^ larger than this. 

Representative finite-source, binary-lens fits are tabu- 
lated in Table 1 a nd shown in Figure El W e recov er the 
fits presented by '.Taroszvhs ki fc Paczvhskil l)2002t) . and 
find many other additional viable fits, with mass ratios 
ranging from g ~ 5 x lO"'^ to ^ 8 x 10"^, at the Ax^ ^ 9 
level. The best binary-lens (and best overall) fit has 
X^ — 130 for 96 dof. If one forces x^/dof = 1 for this 
fit, then mass ratios in the range q = lO"** — 10^^ arc all 
consistent with the data at the ~ 3cr level. 

There two basic classes of fits. In the first class the 
trajectory crosses the planetary caustic, with the deviant 
data point at td occurring when the source is interior to, 
or near, the caustic. The second class of fits pass out- 
side the planetary caustic, with the point at td occurring 
when the source is near the ridge of high-magnification on 
the planet-star axis. The second class actually represents 
a continuous degeneracy in q, d, and a, with mass ratios 
in the range 0.01 ^ <? ^ 0.1 for Ax^ ^ 14. From FigureEl 
wc conclude that while many of the fits would have been 
distinguishable with a few additional data points during 
the perturbation, the family of continuous fits produce 
very similar perturbations. The curves in each of the two 
bottom panels of Figure El deviate from each other by 
< 2% for the majority of the perturbation. Distinguish- 
ing between these models would therefore require rather 
accurate or densely sampled data. T he cause of this 
uncertainty in q is briefly discussed by ijGaudi fc Gouldl 

5. IMPLICATIONS AND DISCUSSION 

OGLE-2002-BLG-055 can be reasonably 
well-fit by large number of different models. 
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iJaroszvnski fc Paczvnskil l)2002() demonstrated that 
the asymmetry exhibited by the primary Hght curve 
can be explained by a simple uniform acceleration, due 
to either parallax deviations arising from the motion 
of the Earth, or non-uniform motion of the source due 
to the presence of, e.g., a binary companion. We have 
demonstrated that the short-timescale deviation can be 
explained by a binary companion to the lens with mass 
ratio 10""* ^ g < 10^^, consistent with a planet, brown 
dwarf, or low mass star. The short-timescale deviation 
can also be explained by a binary companion to the 
source with a flux ratio of ~ 1%. 

Which of these models is most likely? 
IJaroszvnski fc Paczvnskil l)2002j) showed that the 
parallax and constant acceleration fits had essen- 
tially identical x^-vahies for the same number of 
degrees-of-freedom. Thus either are equally viable in 
a goodness-of-fit sense. The inferred timescale of the 
event spans the range = 50— 100 days, which is in the 
range where weak parallax effects are neither surprising, 
nor necessarily expected. Furthermore, there has been 
no estimation of the expected rate of binary-source 
events with detectable asymmetries due to acceleration. 
Therefore it is difficult to argue which origin for the 
observed asymmetry is a priori more likely. In regards 
to the deviant data point, a binary- lens model with 
q ^ few X 10~^ provides the best fit to the data, with 
the binary-source model disfavored at Ax^ — 11 for the 
same number of dof, or slightly more than ^ 3cr level. 
However, this difference in is driven primary by a 
couple of data points just after the perturbation, which 
may be affected by the short-timescale variability seen in 
other parts of the light curve. If one normalizes the error 
bars to force x^/dof — 1 for the best model, then the 
binary-lens model is only favored at the Ax^ — 8 level. 
The inferred blend and source fluxes can also provide 
discrimination between models: one expects these fluxes 
to trace, on average, the distribution of fluxes of un- 
lensed stars in the field near the line-of-sight. The best 
binary-lens models have source and blend magnitudes 
of /s ~ 17 and h ~ 19, whereas the binary-source fits 
yield ~ 18 and If, ~ 17. Unfortunately, without color 
information, it is difficult to definitively distinguish 
between these two scenarios. It is also difficult to argue 
which model is a priori more likely, since the frequency 
of planetary companions to microlens stars is not known, 
and defining the appropriate detection criterion for the 
current event, which was culled by eye from an ensemble 
of alerted microlensing events, is nebulous at best. 

The binary-source model does have the advantage of 
simplicity: a bound companion to the primary source 
can explain both the short-timescale variability and the 
acceleration needed to produce the global asymmetry. 
However, in order for this scenario to hold, the projected 
acceleration vector must be aligned with the projected 
binary-source separation vector. In the particular pa- 
rameterization of the binary-source model used here (see 
this angle AO between the two vectors is essentially 
a free parameter, and thus it can be used to test the 
model. We determine the distribution of A9 from the 
two Markov chains corresponding to the two degenerate 
(in the sense of x^) binary-source models, SL.BS.CA.n 
and SL.BS.CA.p. The second model is immediately ruled 
out, as the primary's acceleration vector is pointing away 



from the companion, with A6 > 90°, wildly inconsistent 
with any bound orbit. The first model has A9 = 36.212^0 
degrees. At first sight, this may also seem very inconsis- 
tent (> 15(t!) with a bound orbit {A9 = 0), however 
there are several reasons why this is not necessarily cor- 
rect. First, although the la uncertainty on A9 is small, 
there is a substantial non-Gaussian tail toward smaller 
values. Enforcing A6 = in the fit yields x^ = 156.7, or 
Ax^ = 15.0 for one less dof. Thus, the bound model is 
only 'ruled out' at the 3.5(t level. Furthermore, as we 
argue below, the period of a bound binary source with 
separation 6 < 1 is likely to be only a few times longer 
than the timescale of the event. Therefore one expects 
significant deviations from uniform acceleration that are 
not accounted for in the simplified model adopted here. 
Thus the inferred value of the A6 may be the result of 
inadequacies in the model. It is therefore conceivable 
OGLE-2002-BLG-055 could be fit by a bound binary- 
source model in which the effects of rotation are included 
self-consistently. 

Accepting the binary-source model as correct, one 
can use the measured parameters of the viable model 
SL.BS.CA.n to determine the Einstein ring r adius pro- 
jected on the source plane te ~ -Dos^^e (Ha n fc Gould! 
fl997l) . where Dos is the distance to the lens. Assum- 
ing circular orbits, the acceleration of the primary due 
to the secondary is simply a — Gm2/cP, where d is the 
semi-major axis of the orbit and m2 is the mass of the 
secondary. The projected acceleration a± is therefore. 



Gm2 



(l-sin^0sin^i)^/^ 



(10) 



where i is the inclination of the orbit, and (p is the phase 
of the primary orbit referenced to the intersection of the 
orbit and the plane of the sky. The constant-acceleration 

fit yields the parameter A = ^ A^ + Ay , which is given 

by, 



A 



(11) 



The binary-source fit yields to, 2- This can be combined 
with the parameters of the primary event to derive the 
two components the projected binary-source separation: 
= 0.5(A/gs) Ar^ - At and by = 0.5(^^/9^) Ar^ - uq. 
Here At = (to, 2 ^ to)/^E, and we have assumed that 
uo 2 1 • Then b can be related to the semi- major axis 

by, 

d = fEfe(cos2 + cos^ i sin^ 0)"^^^. (12) 

Finally, Equations (|10|1 . Hll|) . and H12|l can be combined 
to arrive at an expression for in terms of mostly known 
quantities. 



b^A ^ 



)sin2j)i/2 



(cos^ (j) 



cos^ i SIT? ( 



1/3 



. (13) 

We now use the probability distribution of the binary- 
source model fit parameters, together with equation (fT!I|l . 
to construct the probability distribution for te. Here 
the virtue of MCMC becomes clear: in order to deter- 
mine the probability distribution of te, we simply have 
to determine the value of te at each link in the chain. Be- 
cause the density of points in the binary-source MCMC 
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parameter chain is proportional to the posteriori prob- 
abihty distribution, the resulting distribution of values 
of te is simply proportional to the desired probability 
distribution. 

We first account for the parameters in equation H13|l 
that are not constrained by the binary-source model; 
these include the inclination i, the phase (j), and the mass 
ratio Qs (which enters via b). For each link in the chain, 
we randomly choose ten values of 4> and cosi, assuming 
a uniform distribution of each. We determine Qs using 
the fluxes of the primary and the secondary. We draw a 
source distanc e Dps from the standard Galactic model of 
IHan fc Goull l|l9Ml200l . If the source is in the bulge, 
we assume that the primary is a giant, and therefore has 
a mass of Mi ^ Mq. If the source is in the disk, we de- 
termine the absolute magnitude of the primary from the 
value of Dos J and assuming an extinction of Aj ~ 1.57 
(see 12). We determine the absolute magnitude of the 
secondary from the flux ratio e, primary flux Fs, -Dos, 
and Aj. We estimate the masses of the primary and sec- 
ondary from their absolute magnitudes using an analytic 
approximation to the main-sequence mass-luminosity re- 
lationship derived from the solar metallicity isochrones 
of BcrtcUi eTal] ljl994j) . This then yields the mass ratio 
Qs- We then determine ^e. 

In addition, we also determine the radius -R*^2 of the 
secondary from its absolute magnitude, also usir ig the 
solar-metallicity isochrones of iBertelli et al.l lll994|) . We 
require that yO,_2 = R*,2/rE < U2{td), otherwise we dis- 
card the inferred value of te. We note that the resulting 
distribution of te is not very sensitive to the adopted dis- 
tribution of source distances, for several reasons. First, 
the distance enters through the factor of since b 

formally depends on the mass ratio Qs- However, inspec- 
tion of the expressions for b^ and by reveals that qs mod- 
ifies the acceleration terms, which themselves are oc r|, 
where T2 = (to, 2 — to)/^E- For the binary-source fits pre- 
sented here, t| ^ 10%, and thus these terms are small. 
Second, for bulge sources, where qs is affected the most 
by the assumed source distance, the dispersion in source 
distances is relatively small. Conversely, for disk sources, 
where the dispersion in source distances is large, the in- 
ferred value of Qs is relatively insensitive to the source 
distance (since both primary and secondary are assumed 
to be on the main sequence). Had we simply assumed 
that all sources were in the bulge with Dos — 8 kpc, the 
median inferred value of te would change by < 2%, or 
- O.lcr. 

Figure 13 shows the resulting distribution of te. The 
median and 68% confidence interval is te — 1.87 ± 
0.40 AU. A d optin g the standard Galactic model of 
IHan fc Gouldl l|1995l 1200.'^) . we use the inferred values 
of te and ^e to constrain the mass M, distance Dqi and 
relative proper motion /itrei of the lens. The resulting 
distribution of M (assuming a uniform prior in linear 
mass) is shown in Figure 13 The median and 68% con- 
fidence interval is logM/Mg = -0.28to:42- Similarly, 
we infer Doi = 7.1 ± 0.9 kpc, 6'e = 240;4o/ias, and 

jjL^ei = 5.45tJ;g5 km s~^ kpc~-^. All of these parameters 
are typical for disk-bulge or bulge-bulge lensing events, 
except for /iroi, which is a factor of ~ 6 smaller than 
the median distribution for bulge-bulge events and a fac- 
tor of ~ 8 smaller than the median of disk-bulge events. 
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Fig. 3. — In each panel, the curves show the relative likelihood 
of the given parameter, as inferred from binary-source, constant 
acceleration fit SL.BS.CA.n (see Table 1 and JSj. All curves have 
been normalized so that the likelihood is unity at the peak. The 
solid vertical line shows the median of the distribution, while the 
dotted lines show the 68% confidence limits, (a) The relative like- 
lihood of fgi the Einstein ring radius of the lens projected on the 
source place, (b) The relative likelihood of the mass of the lens, as 
inferred from the measured values of , the tim escal e f e , and the 
standard Galactic model offHan & Gould 1 1995, 200^. (c) The rel- 
ative likelihood of the period of the binary-source in days, (c) The 
relative likelihood of the semi-amplitude K of the radial velocity 
of the primary of the binary source, in kilometers per second. 



This may be due to the fact that the lens and source 
happen to have usually small relative velocities, it may 
indicate that the event is due to disk-disk lensing, or it 
may simply be because the model is incorrect. 

We also infer an Einstein ring radius projected onto 
observer plane of te = IS.Sty'^g^ AU, indicating that, for 
the binary-source model, parallax effects are likely to be 
negligible. 

The probability distributions of several other interest- 
ing quantities can also be constructed from the Markov 
chain. The semi-major axis of the binary is d — 0.80 ± 
0.06 AU. In order to determine the period, we assume 
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a primary mass of mi — Mq, as expected for bulge gi- 
ants. The absolute J-band magnitude of the secondary 
is Ml = 6.60 ± 0.46, which corresponds to a mass of 
7712 = 0.67 ± 0.05 Mq assuming that it is on the main- 
sequence. These corresponds to a late K-dwarf. The 
primary has an absolute magnitude of Mj = 2.01 ± 0.49, 
making it likely to be a A/i ~ Mq bulge giant. The total 
binary mass is M^ = mi + 7712 ~ 1.7Mq. The distribu- 
tion of the binary period is shown in Figure 13 the me- 
dian and 68% confidence interval is P = 200 ± 20 days. 
The binary period is only a factor of ~ 2 larger than 
the microlensing event timescale ^e- This implies that 
the assumption of a uniform acceleration is almost cer- 
tainly violated, and that the binary-source model is not 
internally self-consistent. The constant acceleration ap- 
proximation may not be so bad, as the deviations from 
uniform acceleration will only become significant in the 
tails of the event, when it is not highly magnified. Never- 
theless, the parameters derived from the above analysis 
should be interpreted with caution. 

It is also possible to predict the semi-amplitude of the 
radial velocity of the primary. 



K = 



Gr. 



dM, 



1/2 



(14) 



The distribution of K is shown in Fi 



igureEl 



and 68% confidence interval is IA.O'^'^'q km s 



The median 
^. Radial 

velocity precisions of a few km s^^ should be attainable 
for this source, which has / = 17.9. Thus it may be pos- 
sible to directly confirm the binary-source interpretation 
of this event. 

6. CONCLUSION 

Microlensing event OGLE-2002-BLG-055 exhibits a 
slightly asymmetric, smooth light curve with a sin- 
gle data point that deviates by ^ 0.6 magnitudes 
from neighboring points se parated by several days. 
iJaroszvnski fc Paczvnskil l)2002) concluded that the sim- 
plest interpretation of OGLE-2002-BLG-055 was an 
event with parallax deviations arising from the motion of 
the Earth, and a short-timescale deviation due to a bi- 
nary lens with mass ratio q = 10~^ — 10"^, thereby mak- 
ing OGLE-2002-BLG-055 a candidate planetary event. 

Here we demonstrated that OGLE-2002-BLG-055 can 
be reasonably well-fit by several different classes of mod- 
els, and a wide range of parameters within each model 
class. We found that the data can be fit by many dif- 
ferent binary-lens models whose mass ratios span three 
orders of magnitude, from q = 10"'* to 10~^, thereby 
making the secondary consistent with a planet, brown 
dwarf, or M-dwarf for reasonable primary masses. A 
subset of these binary-lens fits form a family of continu- 
ously degenerate models, whose mass ratios differ by an 
order of magnitude. Astonishingly, the light curves of 



these models differ by < 2% for the majority of their du- 
ration. A binary-source model is also consistent with the 
data, for a secondary /primary flux ratio of ^ 1%. This 
model also naturally explains the global asymmetry of 
the lightcurve as due to the acceleration of the primary 
induced by the secondary. Under the assumption of a 
bound binary-source, this model yields an estimate of 
the Einstein ring radius projected on the source plane of 
fE 1.87 ±0.40 AU. 

All of these fits differ by < Scr, and are essentially 
indistinguishable when the scatter due to likely source 
variability is considered. Unfortunately, the lack of color 
information during the event precludes the discrimina- 
tion of models based on the positions of the source and 
blend on a color-magnitude diagram. 

Although th e primar y go al of the study by 
IJaroszvnski fc P aczvnskil l)2002|) was to affect a modifi- 
cation of the OGLE observation strategy to ensure good 
sampling of short-duration perturbations, rather than ar- 
gue that OGLE-2002-BLG-055 was a bona fide planetary 
event, it is still somewhat disturbing that many binary- 
lens fits were missed, and the possibility of a binary- 
source interpretation was not discussed at all. OGLE- 
2002-BLG-055 serves as an extreme reminder of the de- 
generacies inherent in microlensing events, and high- 
lights the difficulties in interpreting poorly-sampled and 
weakly-perturbed events. These difficulties become es- 
pecially important when attempting to detect planets 
with microlensing. Here observers and modelers need 
to be especially vigilant: in order to produce convincing 
and reliable planetary detections, it is essential not only 
to achieve dense and accurate photometry of planetary 
perturbations, but also to acquire as much auxiliary in- 
formation as possible, and to perform detailed, careful, 
and thorough modeling, in order to ensure that planetary 
detections are robust. 
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Table 1. Model Fit Parameters 



Model 


to 




tE 


h h 


Ax 


Ay 








(HJD-245QQQ.) 




(days) 












SL.SS.CV.l 
SL.SS.CV.3 


2401.67±0.22 
24Q1.61±Q.24 


0.34±0.03 
Q.33±Q.Q3 


68±4 
7Q±4 


17.28 ± 0.12 17.78 ± 0.22 
17.33 ± 0.13 17.70 ± 0.20 








179.2/100 
166.5/98 



SL.SS.CA.l.n 2400.57±0.28 

SL.SS.CA.l.p 2400.57±0.29 

SL.SS.CA.S.n 2400.20±0.33 

SL.SS.CA.3.P 2400.20±0.33 



0.29±0.05 
0.29±0.06 
0.22±0.04 



70^- 17.55^-- 

qO+20 ,7cn+0.31 

17.86t°;t° 
0.22±0.06 1121^^ 17.84t°i8 



± I .'±-J_(j 22 

17 48+''-''i 

,7 nn+0.08 
l(.ZZ_o 

17 24+°-™ 



-0.17 

-0.40 
-0.22 



0.06 
+0.12 
0.20 
+0.05 
0.10 



-0.75 
1.28 
-1.04 



-0 t;R+o-2i 



0.30 

+0.68 

-0.42 
+0.18 
0.57 



1 89+''^'* 

-^■""-O.TO 



148.5/98 
148.5/98 
132.8/96 
132.8/96 



SL.BS.CA.n 2400.00±0.32 0.23±0.05 82 
SL.BS.CA.p 



2400.00±0.35 0.21±0.05 115^^^ 17.97153-^* 



1 7 0(1+0.32 ,7r,,,+0.19 
l/.oy_Q28 ^'-^^-0.60 



17 ,0+0.25 
l'l^-0.55 



1 +0.05 
^-0.06 
-0 ■^■;+0-22 



-0.21: 



0.98^215 0.011+° »^ 2425*0°9+n5 
2425.00lgj« 



, qo+0.96 r. ooq+0.003 
l-^J-0.77 "■'"'^-0.004 



«0,2 

002+°'°°' 



141.7/96 
141.9/96 



















d 




Q 




Q 








BL.SS.CA.a 


2401.17 


0.21 


121.0 


17.95 


17.18 


-0.43 


2.07 


1.19 


4.83 


X 10- 


-4 


0.928 


4.1 X 10- 


-3 


136.8/96 


BL.SS.CA.b 


2401.58 


0.44 


64.9 


16.88 


19.07 


-0.19 


0.59 


1.36 


2.77 


X 10- 


-3 


0.942 


6.6 X 10- 


-3 


133.1/96 


BL.SS.CA.c 


2403.99 


0.38 


56.3 


16.98 


18.51 


-0.09 


-0.54 


1.52 


3.75 


X 10" 


-2 


0.744 


6.3 X 10" 


-3 


130.9/96 


BL.SS.CA.d 


2400.92 


0.41 


54.7 


17.00 


18.43 


-0.12 


-0.44 


1.34 


1.81 


X 10" 


-3 


0.693 


6.3 X 10" 


-3 


133.5/96 


BL.SS.CA.e 


2401.22 


0.33 


82.6 


17.32 


17.71 


-0.26 


0.97 


1.28 


1.51 


X 10- 


-3 


0.943 


5.4 X 10" 


-3 


135.3/96 


BL.SS.CA.f 


2401.24 


0.36 


59.4 


17.20 


17.90 


-0.13 


-0.52 


1.30 


2.03 X 10- 


-3 


0.703 


5.7 X 10" 


-3 


135.1/96 



Note. — SL=smgle lens, BL=bmary lens, SS=slngle source, BS=bmary source, CV=constant velocity, CA— constant acceleration, l=one data point 
removed, 3— three data points removed, n— negative constant acceleration fit, p— positive constant acceleration fit, (a, b,c,d,e,f)— binary-lens models shown in 
Figure 2. 
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